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Abstract 

The shape of the electron energy spectrum in 3 H /3-decay permits a direct assay of the absolute 
scale of the neutrino mass; a highly accurate theoretical description of the electron energy spectrum 
is necessary to the empirical task. We update Sirlin's calculation of the outer radiative correction 
to nuclear /3-decay to take into account the non-zero energy resolution of the electron detector. 
In previous 3 H /3-decay studies the outer radiative corrections were neglected all together; only 
Coulomb corrections to the spectrum were included. This neglect artificially pushes rr? v < in a 
potentially significant way. We present a computation of the theoretical spectrum appropriate to 
the extraction of the neutrino mass in the sub-eV regime. 
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Empirical evidence of neutrino oscillations in atmospheric, solar, and reactor neutrino data ^[^UHl 
compels the existence of non-zero neutrino masses, yet such experiments are insensitive to the absolute 
scale of a neutrino mass, for the oscillation experiments determine Am| = mf — m 2 , where nii is the 
mass of neutrino i. To determine the absolute value of the neutrino mass requires different methods. 
Cosmological constraints on the neutrino mass do exist [I], though our focus shall be on the study 
of the electron energy spectrum in tritium /3-decay near its endpoint, as this represents the most 
sensitive terrestrial measurement. The spectrum shape constrains the mass of the neutrino, be it of 
Dirac or Major ana character, and the inferred mass is insensitive to phases in the neutrino mixing 
matrix — in contradistinction to the constraint on the neutrino mass from neutrinoless double /?- 
decay. An accurate theoretical description of the expected electron energy spectrum is crucial to the 
determination of the neutrino mass; this demand grows as the sensitivity of the experiments increase. 
Indeed, future studies expect to probe the neutrino mass at the sub-eV level 0- It is our purpose to 
realize a theoretical form of the requisite accuracy, though we shall begin by describing the form used 
in earlier tritium experiments. 

With an anti-electron neutrino of mass m v , neglecting neutrino mixing for simplicity, the Fermi 
form of the electron energy spectrum for tritium /3-decay is jH] 
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where Gf is the Fermi constant, p e , E e , and .E™ ax are the momentum, energy, and maximum endpoint 
energy, respectively, of the electron, and |A4| 2 is the absolute square of the nuclear matrix element, 
with \Ai\ 2 ~ 5.3. A form of this ilk has been used to bound m u in previous experimental analyses 
of molecular tritium /3-decay [7J El El HB EEU El EC3 El- Following the usual practice, we include a 
non-zero neutrino mass in the phase space contribution only. We set h = c = 1 throughout. The 
Fermi function, F(Z, R e , E e ), captures the correction due to the Coulomb interactions of the electron 
with the charge Ze of the daughter nucleus ^] • We adopt the usual expression JHj > derived from the 
solutions of the Dirac equation for the point-nucleus potential —Za/r evaluated at the nuclear radius 
Re C3i ^ differs from unity by a contribution of 0{a). The Fermi function includes the dominant 
electromagnetic effect, though an accurate extraction, or bound, of the neutrino mass does demand 
the inclusion of the remaining 0(a) correction. We shall demonstrate this point explicitly. This last 
effect, termed the radiative correction, is conventionally separated into an "inner" piece which is 
absorbed in \M.\ 2 , as it is energy independent and thus of no consequence to our current study, and an 
"outer" piece 5r [THj . The outer radiative correction applied in /3-decay studies is due to Sirlin 
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with M and m e the proton and electron mass, respectively, and L(x) the Spence function. As in 
Ref. JH]) we neglect terms of relative order a(E e /M)ln(M/E e ), aE e /M, and smaller, throughout. 
The total 0(a) correction is given by F(Z,R c ,E e ) - 1 + (a/2ir)g s (E e , Ef 3 *). Note that S R results 
from averaging {a/2Tr)gs (E e , E^ ax ) over the electron energy spectrum. To assess the relative sizes 
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of F and gs, we note that in the endpoint region of tritium /3-decay, for E e — m e = 18.5 keV, e.g., 
F ~ 1.19, whereas (a/2ir)g s ~ 0.02. 

Since the absolute neutrino mass scale is inferred through the shape of the electron energy spectrum 
in the endpoint region, it is crucial to predict the shape of the theoretical spectrum with high accuracy. 
To this end, we update the calculation of Ref. ^H] to take the energy resolution of the electron detector 
into account. To understand the significance of this, we recall that Sirlin's function contains not only 
virtual photon corrections to the /3-decay process but also bremsstrahlung contributions, to yield an 
additional real photon in the final state. Only their sum is infrared finite; the infrared divergence in 
each contribution is regulated by giving the photon a small mass A, with the A — > limit to be taken 
after the sum has been computed and the infrared divergent pieces cancelled. The finite portion of 
the bremsstrahlung contribution is sensitive to the precise manner in which the experiment is effected. 
In Sirlin's function the energy resolution of the electron detector is implicitly assumed to be zero; 
that is, the e~ and 7 are always distinguishable. A consequence of this is that Eq. (|HJ) contains a 
logarithmic divergence as E e — > E™ ax . This singularity can be removed by including soft photon 
contributions to all orders in perturbation theory |19j : however, it can also be removed by taking the 
energy resolution of the detector into account, as we consider here. The bremsstrahlung contribution 
comes from integrating the photon energy over its entire kinematic range, namely 

^\dE-) 9S > b ^ E e ] -J x H dE e dE, J' (4) 

where d 2 Y 1 {E e ,E 1 )/dE e dE 1 is the doubly differential decay rate for the radiative /3-decay process. 
Nevertheless, the detector energy resolution AE can, in principle, influence the shape of the electron 
energy spectrum. That is, for E 7 < AE, the electron and photon cannot be distinguished; indeed, this 
is precisely why the bremsstrahlung contribution can enter to render an infrared- finite, 0(a) radiative 
correction to /?-decay. In this event the detector records the sum of the electron and photon energies 
as the "electron" energy. For E 7 > AE, however, the electron and photon energies are distinguishable, 
and their energies can be recorded separately. To separate the contributions we note that the total 
bremsstrahlung contribution to 5r must be insensitive to such experimental details. In passing, we 
note related discussions of the impact of the detection threshold for bremsstrahlung photons in the 
radiative corrections to v capture on deuterium [201 Ell 1221 G3] , as well as to v — e scattering [21] . The 
bremsstrahlung contribution to the total outer radiative correction 5r — we retain the photon mass 
A throughout — is 

S R ,b= I dE e j E dE 7 /(E e ,E 7 ), (5) 

where f(E e , E 7 ) = <i 2 r 7 (E e , E 1 )/dE e dE 1 . To find 5 R we must divide <5r & by a normalization factor N , 
determined by integrating Eq. Q over the allowed phase space with F = 1 and m v = 0. However, if the 
detector energy resolution is "infinite," that is, if the electron detector always records £■ = E e +E^ |20| . 
the total bremstrahlung contribution can be rewritten as 

/■E e max rS-m e 

hb = dS dE^ f{£ - E y , . (6) 

Jm e +X J\ 

In both Eqs. © and Q, the integration over EL yields the bremstrahlung correction to the electron 
energy spectrum. Although 5R t b is universal, the shape correction to the electron energy spectrum 
is not. Let us now determine the shape correction for a finite energy resolution AE. That is, for 
E 7 < AE, £ = E e + E 7 is recorded, whereas for E 7 > AE, E e is recorded. We can reorganize the 
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total bremsstrahlung contribution in the following way: 
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Note that letting AE -> A yields Eq. ©, whereas letting AE -> E™ 3 * - m e yields Eq. ©. If AE 
has a non-infinitesimal value, then infrared divergences are restricted to the first two terms — in the 
third, we may let E 1 = \/k 2 + A 2 — ► A; with impunity. Thus we introduce 



a dT 



and 



to realize 



2?r d£ 

a dTo 
2tt1eZ 



..max 

J A (A: max , £) = y d£ 7 f{£ - E y , 

l(AE,E e )= / E " dE 7 f(E e ,E 7 )^ / £ ^dkf(E e ,k) (8) 

jAJ? JA£ 



<? 6 (A£, E e ,E?°*) = Q(AE + m e - E e )l x (k™ x = E e - m e , E e ) 
+ G(E e - (AE + m e ))T A (fc max = AE, E e ) 
+ G(£ e max - AE - E e ))l(AE,E e ) . (9) 

With this, we determine that the radiative correction to the electron energy spectrum is 

g(AE, E e , Er x ) = g b (^E, E e ,E^) + g v (E e ) , (10) 

where the virtual photon contribution g v = 2 A — 3/4, with A as reported in Eq. (16) of Ref. 23 . To 
compute the integrals of Eq. © and thus gb(AE, E e , E™ ax ) we adapt the computation of radiative 
neutron /3-decay in Ref. [251 to this case. In specific, we use the absolute squared matrix elements of 
Eqs. (13-16), and Eq. (20), dividing the latter by 8M 2 , in that work to determine f(E e ,E 7 ). In doing 
the integrals, we can neglect the recoil corrections to the phase space, so that E u = £™ ax — E e — k. 
As a result, the dependence on the weak, hadron coupling constants is captured by g 2 ? + 3g\, which 
we replace by |-M| 2 ; we also update the mass factors which appear as appropriate. We have verified 
that our computation of gb(AE = 0, E e , E™ ax ) using this procedure and Eq. (JSJ), in concert with g v of 
Ref. [23!) yields Shim's result JE]) Eq. (j3J. For finite AE we find the following results: 
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where j3k = yJS — k) 2 — m 2 / (£ — k). (We note that I£ can be brought to closed form via the 
substitution t = though we omit the resulting expressions here.) Moreover, 



1(AE, E, 

m 



-4 In 



El 



E„ 



+4 



max 
e 



E, 



AE 

(AE) 2 ) 



max 



E, 



AS) , AE- 



E„ 



+ 



E e - AE) 



max 



E, 



,((e, 



max 

e 



E P 



(AEf 



E e (E, 



max 
e 



+-tanh _1 /3 

P 



2 In 



E e ) 



E e 



AE 



+ 2 



(E, 
(E, 



max 

e 



-E e ) 2 
E e - AE) 



{El 



E e — AE) A(E» 



E P 



E e 

(AE) 2 ) 



+ 



3 

m 



max 

e 



E e 



(AEf 



E e (Ef™-E e 



Ee 



(AEf 



(er - e c ) 

((Ef^ - Ee) 2 - (AEf) 



(E* 



E P . 



2 ((£ e max - 
3 E e (E, 



im 

-E e ) 3 
E e f - (AEf 



2E 2 
E e f - (AEf ) 



El(El 



Eef 



1 ((B e max ~ 
4 E 2 (E, 



-E e ) 
E e f 



(AEf 



Eef 



(13) 



We note that the In A term in I\(k ma,x , £) cancels the concommitant infrared divergent term in g v to 
yield a finite result in the A — > limit, irrespective of the detector energy resolution AE. Using these 



formulae, we find our g(AE = E T e ' 
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is in accord with the result of Vogel, Ref |2f)j . 



We report g(AE, E e , £'™ ax ) in the endpoint region of tritium /3-decay, which results from Eqs. (jSJ- 
H3J), in Fig. HI We have verified that the integration of (dT /dE e )g(AE, E e , E^ 11 ^) over E e yields a 
universal value of 5jib for all AE. The inclusion of a finite value of AE removes the logarithmic 
divergence in Sirlin's function as E e — ► -Eg iax . The numerical shifts associated with the inclusion of 
the AE dependence generally are crudely comparable in size to the leading (D(Za 2 ) correction |2fi| . 
— Za 2 m(M/m e ), though the latter, of course, contains no E e dependence. 

We can now proceed to evaluate the changes these theoretical corrections make to the shape 
assumed in previous experimental assays of the v mass. For definiteness, we also include the leading- 
order recoil corrections to the electron energy spectrum. We may once again adapt the results from 
neutron decay to this case. We adopt the notation of Bender et al. in Ref. j^Tj, though the couplings of 
the hadronic weak current are now nuclear form factors evaluated at zero momentum transfer. Noting 
Ref. 122], we replace the absolute, squared nuclear transition matrix element, which we have taken to 
be |^| 2 = |A^ | 2 = ^ + 3<7i, with 
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Note that M A = 2809.4319 MeV EH1 is the tritium mass and £f ax - m e = 18.57 keV 14 , so that 
^max jy[ A ^ i Q . io -4 . Since |.Mo| 2 can be absorbed into the overall normalization of the decay rate, 
the function 1Z represents the first appearance of nuclear-structure effects in the prediction of the 
electron-energy spectrum in tritium /3-decay. The form factors which enter are largely determined by 
the symmetries of the Standard Model (SM), so that the subsequent uncertainty in the predicted recoil 
correction, which is itself of small numerical size, is very small. In writing Eq. (|15|) . we have assumed 
the validity of the conserved-vector-current (CVC) hypothesis and have neglected the form factors 
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Figure 1: The "outer" radiative correction g(AE, E e , i?™ ax ) as a function of the electron detector 
resolution AE and the electron energy E e in the endpoint region of tritium /3-decay. The solid line 
is Sirlin's result ^Hj, for which AE = 0. The dotted line corresponds to Vogel's result [20], for 
which AE = E™ ax — m e . The remaining curves correspond to AE = 1 6 100, and 1000 
eV, respectively, moving in sequence from the solid line to the dotted one. The chosen values of AE 
correspond to those of the planned and recent experiments to which we refer. 



associated with second-class-current contributions. In the context of the SM, this is tantamount to 
neglecting the effects of isospin violation, so that the recoil term is subject to corrections of 0(1%). 
The vector coupling gy is also unity by the CVC hypothesis; the computed correction due to charge- 
symmetry breaking in the overlap of the 3 H- 3 He wave functions, due to Towner, is 5 C = 0.06% |3Uj : 
we note g v = g'y 2 (l — S c ), where g' v 2 absorbs the inner radiative correction A R and |14d| 2 , with V u d 
a CKM matrix element. The CVC hypothesis determines the weak-magnetism coupling fi from the 
measured 3 H and 3 He magnetic moments j2Sj, to yield jijg'y = —3.0533; we ignore the possibility of 
a inner radiative correction idiosyncratic to fi- The 3 H half-life determines g v + 3g^ up to corrections 
of recoil order; this is sufficient to determine the couplings which appear in the recoil-order expression. 
In specific, we have M (1 - *c + ^\ I d'y 2 )' 1 = G 2 F g' v 2 f{\ + S R )t 1/2 /K with g' y 2 = (1 + A R )\V ud \ 2 
and K = 2vr 3 ln2/m^. We use G F , a, m e , and h as given in Ref. [HI], A R = 0.0240, V ud = 0.9740 as 
in Ref. [HU, an d the half-life ti/2 = 12.3 yrs as recommended in Ref. 28 j. Finally we use the integral 
of F + (a/2Tr)gs, noting Eq. (4) of Ref. ^Hl for F, over the allowed phase space to fix /(l + 5r), for 
which we find 2.9109 • 10~ 6 . We use R e = 1.68 fm throughout [HHH]- This yields g A /g'y = 1-22; note 
that this ratio of couplings implicitly contains the quenching of the Gamow- Teller matrix element due 
to nuclear structure effects. 

Armed with these results, we now proceed to evaluate the change in the electron energy spectrum 
upon the inclusion of the outer radiative correction, g(AE, E e , i?™ ax ). We illustrate A(dT/dE e ) = 
{cLTq/ dE e )ag{AE, E e , £^™ ax )/(27r) in Fig. [21 using the energy resolution of the Mainz experiment |14j . 
AE = 4.4 eV. In this figure, the inclusion of the AE dependence is of little impact; the resulting 
curve is hardly distinguishable from that which results from the use of Sirlin's function, Eq. ©. 
The recoil corrections are included as well, so that we employ dT/dE e = dT^/ dE e (l + 1Z + F + 
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Figure 2: The change in the electron energy spectrum, dT/dE e , upon the inclusion of g (AE, E e , E™ 3 *-), 
as a function of E e . The solid line has AE = 4.4 eV as per Ref. |14j : the dashed line also includes recoil 
corrections as per Eq. (|14|). For reference the change in the theoretical form of the energy spectrum 
assumed in Ref. |14j . if m. 2 = 0^—4 eV 2 , is shown as the dot-dashed line as well. 



aj {2i:)g(AE, E e , E™ ax )); they are rather small, though they are appreciable. The analysis of Ref. 
assumes the theoretical form given in Eq. inferring m 2 = —3.7 ± 5.3 ±2.1 eV 2 from their exper- 
imental data. Thus, for reference, we also show A(dT/dE e ) = dT F (ml = -AeV 2 )/dE e - dT F {ml = 
OeV 2 ) / dE e in Fig. [21 Note that employing m 2 > yields a A(dT / dE e ) which differs in sign from that 
generated with g(AE , E e , E™ ax ) . It is apparent that the change in the theoretical energy spectrum 
due to the neglected 0(a) correction acts to increase the electron energy spectrum; this effect is also 
realized through a negative value of m 2 in Eq. (^Q). The neglect of the outer radiative correction 
generates a negative shift in to 2 . We emphasize that this shift is a consequence of the change in 
shape in the electron energy spectrum induced by the outer radiative correction. We illustrate this 
in Fig. 13 in which we show the ratio of the corrected to uncorrected Kurie plots, recalling that the 
Kurie plot is K(E e ) = [(dT /dE e )(l/Fp e E e )} 1 / 2 versus E e . This ratio is simply (1 + g eE ) 1 / 2 to O(o), 
where g eS = K + (a/2ir)g{AE, E e ,Ef ax ). 

To estimate the impact of the remaining 0(a) correction on the value of the neutrino mass, we 
introduce the fit function used in earlier work |14j : 

rm e ^- = Am- 4 F(Z, R e , E e ) Pe E e (E™ x - E e )y/ (E™* - E e )* - m* + B , (16) 
arj e 

where A, E™ ax , to 2 , and B are all fit to the electron energy spectrum. For our comparison, which 
we effect for purposes of illustration, we have neglected the contributions associated with the excited 
final states of the daughter 3 He + - T molecule and include the elastic contribution only. Note that 
B represents a constant experimental background, so that we set B = 0. Fitting the remaining 
parameters in Eq. (|16|) to the dashed curve in Fig. [21 yields the comparisons shown in Fig. [1J We fit 
the last 70 eV of the electron energy spectrum and set AE=AA eV, in an attempt to simulate the 
conditions in Sec. 6.3 of Ref. for which to 2 = — 3.7±5.3±2.1 eV 2 was inferred. We fit the quantity 
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Figure 3: The ratio of the corrected to uncorrected Kurie plots, namely (1 + g c s) ' , with g c R = 
TZ + (a/2ir)g(AE, E e , E™ x ), as a function of E e . The solid line has AE = 4.4 eV, the dashed line has 
AE = 1 eV, and the dotted line has AE = 0.1 eV. 



Tm e dF/dE e , where the lifetime r is r = i 1 / 2 /ln2, noting, for reference, that A = 3.43537 • 10 5 in our 
theoretical curves. Using the MIGRAD minimization program in the CERN package ROOT 3*5], we 
find two different fits of comparable quality, which possess very different values of m 2 ; apparently a 
significant shift in m 2 can be accommodated by normalizations A which differ by ~ 2%. In both 
cases m 2 < 0; we cannot successfully fit our curves using a non-negative value of m 2 in Eq. (fTB|) . 
In constrast, if we fit Eq. (J16|) to a curve containing the Fermi function only, these features do not 
occur. In this case, we find fits with |m 2 | < 1 eV 2 are consistent with the input curve; m 2 need not 
be negative definite. 

The refinements we have introduced, namely, the AE dependence to the outer radiative corrections 
cum the recoil corrections, shift ml at no larger than the 0(1 eV 2 ) level. As Fig. |31 makes clear, this 
conclusion is sensitive to the precise value of AE, as well as the interval in E e over which the neutrino 
mass is fit. For other choices of these parameters, their relative impact could be more significant. 
Given the results shown in Fig. @, we cannot make a robust conclusion concerning the absolute scale 
of the shift in m 2 our previously neglected theoretical corrections would induce in a more realistic 
analysis; nevertheless, we can say that the consequence of neglecting these terms is to push m 2 < in 
an artificial way. We presume that with the replacement of Eq. (|16|) with a fit function incorporating 
the radiative and recoil corrections we have calculated such artificial shifts would disappear. In the 
analysis we have effected, this turns out to be the case. 

In this letter we have evaluated the O(a) outer radiative correction, g(AE, E e , .E™ ax ), to the 
electron energy spectrum in 3 H /3-decay, so that {a/2ir)g + F — 1, where F is the Fermi function, 
constitutes the complete 0{a) correction to the electron energy spectrum. We have updated the 
calculation of Sirlin 18 to include the dependence of the outer radiative correction on the detector 
energy resolution AE; as a consequence the 0(a) correction we compute to the shape of the electron 
energy spectrum is finite as E e — ► i?™ ax . However, as necessary, it has no impact on the radiative 
correction to the total decay rate. Interestingly, the outer radiative correction was omitted all together 
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Figure 4: The electron energy spectrum for E™* — E e < 70 eV, as a function of E e . The solid curve 
shows the theoretical spectrum to be fit, which includes g(AE, E e , Eg 3 "*) with AE=4A eV and recoil 
corrections as per Eqs. (|14I15|) . The dashed curve is realized from Eq. (|16|) using _E™ ax = 18.57 keV, 
A = 3.55020 • 10 5 , and m 2 = -67 eV 2 . In constrast, the dot-dashed curve has A = 3.49446 • 10 5 and 



ml = -0.01 eV 2 . 



in earlier studies of tritium /3-decay [7J El El EE EEU El CS1 El E] ; we have shown that the shape 
correction associated with this ~ 2% shift mimicks a negative value of m 2 . We believe it is necessary 
to update earlier experimental analyses to take this theoretical correction into account, to realize an 
accurate determination of the neutrino mass. A highly accurate theoretical spectrum can be found 
by modifying Eq. the form used in earlier experimental analyses of 3 H /?-decay, through the 
substitution F -> F* + (a/2ir)g(AE, E e , E™ ax ) + K, using Eqs. PT3|) and Eq. (JTSJ). Our focus 
has been on g(AE, E e , E™ ax ) and 1Z; theoretical corrections to these terms accrue from i) 0(Za 2 ) 
corrections, which are known and ii) 0(1%) corrections to the recoil-order term, Eq. ()15|) . but 
such corrections would appear beyond the scope of current and planned experiments. The corrected 
Fermi function F* , which includes corrections such as those due to the finite nuclear size and to charge 
screening of the nuclear charge by atomic electrons, is detailed in Ref. [HI]; recoil corrections [221 an< ^ 
outer radiative corrections, as calculated by Sirlin ^H], are considered in this reference as well. The 
outer radiative corrections are the largest of these corrections |34j . Realistic experimental conditions 
demand that Eq. (|16|) . as well as the fit form we advocate, be adapted to include the population of all 
the excited final states i of the daughter He + -T molecule; the excitation energies and amplitudes are 
computed from atomic theory. As a consequence, the elucidation of the neutrino mass in 3 H /?-decay 
relies on atomic physics input which cannot be wholly subjected to exhaustive, independent empirical 
test. Nevertheless, from the viewpoint of the theoretical radiative and recoil corrections, a sub-eV 
determination of the neutrino mass should be possible. 
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